The infrared conductivity of graphene 
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We study the infrared conductivity of graphene at finite chemical potential and temperature taking into ac- 
count the effect of phonons and disorder due to charged impurities and unitary scatterers. The screening of the 
long-range Coulomb potential is treated using the random phase approximation coupled to the coherent poten- 
tial approximation. The effect of the electron-phonon coupling is studied in second-order perturbation theory. 
The theory has essentially one free parameter, namely, the number of charge impurities per carbon, nf. We 
find an anomalous enhancement of the conductivity in a frequency region that is blocked by Pauli exclusion and 
an impurity broadening of the conductivity threshold. We also find that phonons induce Stokes and anti-Stokes 
lines that produce an excess conductivity, when compared to the far infrared value of <7rj = (tt /2)e 2 jh. 
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PACS numbers: 81.05.Uw, 73.25.+i, 72.80.-r 

Although there has been enormous experimental-!- and 
theoretical? progress in understanding the physical proper- 
ties of graphene since its isolation in 2004 3 , the important is- 
sue of graphene transport remains unsettled. More than 20 
years ago, E. Fradkini showed that the D.C. (static) conduc- 
tivity of graphene, at the charge neutrality point, ctd.c. (a 4 = 
0) = a(u) = 0,/i = 0) (oj is the frequency and fj, is the 
chemical potential measured relative to the Dirac point), can- 
not be described within the standard Boltzmann approach of 
metals, because the Dirac-like electronic excitations have in- 
finite Compton wavelength (which is cut-off by the size of the 
sample), violating the assumptions for the validity of Boltz- 
mann transport'. Fradkin showed that the proper way to com- 
pute the conductivity is through the Kubo formula treating 
the impurities in a self-consistent way. The Kubo formula 
predicts a universal, impurity independent, D.C. conductivity: 
fD.C.theo. (a* = 0) = (4/tr)e 2 /h. Nevertheless, experiments 1 
find that crrj.c.,exp.O u = 0) ~ 4e 2 /h with sample-to-sample 
variations by a factor of 2, which importantly are in the direc- 
tion of higher conductivity, i.e., further away from the theo- 
retically predicted value (the so-called mystery of the missing 
7t). This result has been assigned to the macroscopic inhomo- 
geneity and non-local transport in graphene samples 1 . Never- 
theless, there is still no consensus in the theoretical commu- 
nity on the origin of this effect?. 

In order to settle the issue of transport and scattering mech- 
anisms in graphene, an aspect of major scientific and techno- 
logical significance, it is important to study electronic trans- 
port away from the static regime. In this regard, the frequency 
dependent (A.C.) conductivity, er(w,/x) (we use units such 
that H = 1), provides important information on the scatter- 
ing mechanisms of the charge carriers for frequencies ui<2fi. 
The basic physical processes involved in the A.C. conductiv- 
ity are easy to understand. A graphene sample is illuminated 
with light of frequency ui and vanishing small wavevector that 
causes creation of particle-hole pairs (pair creation), as shown 
in Fig. Q] At zero temperature and in the absence of disor- 
der or phonons only particle-hole pairs with energy greater 
than 2/i are allowed since all the states with energy between 
— n and +fj, are forbidden transitions due to Pauli's exclusion 
principle. In this case the A.C. conductivity is simply a step 



function <t(uj, fi) = cto6(cj — 2/i), where <jq = (7r/2)e 2 //i is 
the optical conductivity that has been measured recently^. The 
far infrared conductivity is insensitive to phonons, impurities 
(since these affect only the low energy part of the spectrum) 
and band structure effects when the frequency of incident light 
is much larger than 2/i and much smaller than the electronic 
bandwidth, W (« 9 eV) 7 -. Nevertheless, as we are going to 
show, the infrared spectrum is very sensitive to phonons and 
impurities and the response of the system deviates substan- 
tially from the non-interacting clean problem. 

Phonons of frequency uiq can be either absorbed or emit- 
ted by the Dirac electrons. When these phonons are at the 
center of the Brillouin zone (T point) they can be probed by 
Raman spectroscopy^, playing an analogous role as light in an 
A.C. conductivity experiment, that is, creation of particle-hole 
pairs^. Furthermore, impurities play a fundamental role at low 
energies since it is known that they produce strong broaden- 
ing of the line-widths 9 . We stress once more, as in the case of 
D.C. transport, that the impurity broadening has to be calcu- 
lated self-consistently. 

In this paper we compute a(u!,/j,) taking into account the 
combined effect of impurities and phonons. We assume that 
there is a density nf of charge impurities per carbon which 
might be trapped in the substrate, on top of graphene, or in 
the interface of graphene and the substrate. We model the 
screening of charge impurities via the random phase approx- 
imation (RPA) together with the coherent potential approxi- 
mation (CPA), which gives us the self-consistent density of 
states. We also assume a density m of unitary scatterers that 
exist due to structural disorder (edge defects, cracks, vacan- 
cies, etc). The effect of unitary scatterers is only important 
in producing a finite density of states at the Dirac point and 
this can be obtained with arbitrarily small values of n,. We 
assume throughout the paper that impurities are dilute and the 
structural disorder is very weak, that is, 1 ^> nf n, — > 0. 

We have checked that the effect of in-plane acoustic 
phonons is negligible and they will be ignored in what fol- 
lows. We assume that the coupling of graphene to the sub- 
strate is strong enough to shift the flexural phonon frequencies 
away from the infrared regime, allowing us to ignore them for 
the moment beingiS. Hence, we have kept only the optical 
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The phonon Hamiltonian has the for m 11 ' 12 ! 13 ' 14 ' 15 : 

H ph . = J2 



R 



Pi(R) , P|(R + 5 3 ) 



2AI, 



c 



2M C 



^^[(u A (R)-u B (R + S)).SY 



Particle 




FIG. 1 : (color online) A photon of frequency u creates a particle- 
hole pair around the Dirac point. The electronic states with energy 
between — fj, and +fi are blocked transitions due to the Pauli exclu- 
sion principle. In the figure the phonon frequency uio is assumed to 
be smaller than /i. Light with with energy 2 x (ujq + /i) lead to 
the generation of an anti-Stokes emission. The Stokes emission with 
energy 2 x \fj, — u>o \ lies in the forbidden region. 



phonon modes. The phonon frequency, and the value of the 
electron-phonon coupling is fixed from Raman experiments^ 
and therefore they are not fitting parameters here. In fact, we 
have only one fitting parameter, namely, nf. 

The Hamiltonian has the form: 



H = H n 



H 



P h. 



H , 



- P h. 



imp. j 



(1) 



where 



Ho = -t (4Wa (R + S)+ h.< 



(2) 



R., 



is the nearest-neighbor tight-binding kinetic energy where 
a^(R) 5)) creates an electron on site Rof sub-lattice 

A(B) with spin a (a =|, J,), i(«3 eV) is the hopping energy 
and 6 are the nearest-neighbor vectors 2 -. 



+ E IT [cos(0(H, 5)) - cos(# )] 2 1 , (3) 

where ua,b are the displacements of the A (B) atoms from 
equilibrium {Pa.b the momentum operator), Ale (= 12 a.u.) 
is the carbon mass, a (sa 500 N/m) is the stretching elastic 
constant, 0(R, d) = (9^ is the angle formed between the i—j 
bond and the i — k bond (6q = 120° is the equilibrium angle) 
and (3 (~ lON/m) is the in-plane bending elastic constant 
(a = 1.42 A is the carbon-carbon distance). Although (0) 
describes both acoustic and optical phonon modes, we focus 
on the optical modes which can be written as: 



»(R) = (u A (R) - u B (R + St 
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with frequency lu% = 3{a + 9f3/2)/M c « 0.2 eV (w 1600 
The electron-phonon Hamiltonian can be written as: 



1 dt 



P h 
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B Q ^)[e tk - s aUk + Q)b a (k) 



e- zk - s bi(k)a a (k-Q)} 



(5) 



where dt/da ~ 6.4 eV A is the electron-phonon coupling, 
Bq v creates a phonon of momentum Q, polarization v (po- 
larization vector e,y(Q)), and frequency lo v (Q) (N c is the 
number of unit cells)^ 2 -. 

The impurity Hamiltonian has the form: 

H lmp . = jfY, V{q)[al{p)a a {p + q) + bl{p)b a {p + q)}. 

(6) 

For a screened impurity of bare charge Ze the potential V(q) 
is given 16 by 



V(q) = 



Ze 2 e- qd 
2eA c q + 7 



(7) 



where e = 3.9 is the SiC>2 relative permittivity, d is the 
distance of the impurity to the graphene plane, and 7 = 
p(/x)e 2 / (2e^4 c ) is the RPA screening wavevector- where p(n) 
is the self-consistent density of states (A c = 3v / 3a 2 /2 is the 
area of the unit cell). Unitary scatterers are modeled using a 
local potential V(q) = U and taking U — > 00. 

The effect of a dilute concentration of unitary scatterers can 
be calculated exactly using the T-matrix, leading to a retarded 
impurity self-energy of the forrn^: 



Z k G (k,iu + i0+)/N c 



(8) 
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where Go(k,u>) is the free electron Green's function 
associated with Hamiltonian ©. Since po(u) = 
~ I/ 71 " Efc 3Go(fe, uj)/N c oc |w| is the bare density of states 
of the clean problem, it is easy to see that (O leads to a diver- 
gence of Ejj 11 *'^ — * 0) at the Dirac point which is unphys- 
ical. Hence, the problem has to be treated self-consistently 
by replacing G {k,u)) by G unit -(fe, w) = G (k,w)/(1 - 
Go(fc, w)S^ nlt (w)) in ([8j. In this case, one can show that 
the self-energy becomes finite at the Dirac point 1 -, as shown 
in Fig. [2] leading to a finite density of states at zero energy. 

The self-energy due to charged impurities is calculated in 
second order perturbation theory as: 

Q 

E c (fe, iu n ) = ^ E - P)| 2G V . (9) 
c p 

where a term of the form nf V(0) was absorbed in the defini- 
tion of the chemical potential. The self-energy (0 is depen- 
dent both on the momentum k and on the frequency. However, 
we are interested on the effect of the self-energy for momenta 
close to the Dirac point (q = K = ^(1/3, V3/9)). Within 
this approximation, the imaginary part of the retarded self- 
energy becomes diagonal and momentum independent, read- 
ing (d ~ 0): 

Notice that the imaginary part of the self-energy behaves like 
\ui\ at low frequencies and vanishes as l/\u>\ at large frequen- 
cies (see Fig|2|i. 

Hence, the electron Green's function in the presence 
of impurities is written as: G~ 1 (k,ui) = GQ 1 (k,oj) — 
S unit (w) - S c (K,cj). Notice that the density of states, 
p(u>) = —1/(ttN c ) J2k ^G(k, uj), should be computed self- 
consistently since the screening wavevector 7 in (0 and ( [Tol l 
depends on p(p). 

The self-energy due to electron-phonon interaction is also 
computed at the Dirac point in second order perturbation the- 
ory: 

1 ' " 7 2 \da) M c u N c ^ 
x \ D °(Q> w m)G(K - Q, iu n - iv m ) (11) 

m 

where D°(Q,iv m ) = 2uiq / ((iv m ) 2 — [ujq] 2 ) is the phonon 
Green's function. Notice that G(k,iuj n ) is the impurity 
dressed electronic Green's function. Due to the exclusion 
principle, the imaginary part of the electron-phonon self- 
energy vanishes when p, — ujq < uj < p + uj^, at T = 0. At 
high frequencies the self-energy follows the electronic density 
of states and is, therefore, linear in uj, as shown in Fig|2] 

In the presence of an electromagnetic field the hopping en- 
ergy changes to: 

t -> te %eA ^- s . (12) 
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FIG. 2: (color online) Imaginary part of the electronic self-energy 
due to unitary scatterers, Coulomb impurities, and phonons. The 
impurity concentration is = nf — 2 x 1CP 4 . We have set p = 
OeV (solid line), p = 0.29eV (dashed line) and T = 45 K. 

Expanding the exponential up to second order in the vector 
potential A(t) and assuming the electric field to be oriented 
along the x direction, the current operator is obtained from 
jx = -8H/dA x (t) leading to j x = g + A x (t)j£. The 
Kubo formula for the conductivity is given by: 

.H = ^-(-^[0-f) + A^ + l + )], (13) 

with A s = N C A C the area of the sample and 

A x *(i"n)= [ dTe^(T T £ (t)£ (0)) , (14) 
Jo 

is the current-current correlation function. The finite fre- 
quency part of the real part of the conductivity is given by: 

&H = - / — e(w , ,w)[/(fiw , -/i)-/(fia; , +fiw-A*)], 

7TAI J UJ 

(15) 

where f(x) is the Fermi function and Q(u>', uj) is a dimension- 
less function that depends on the full self-energy and will be 
given elsewhere^. The main features of the conductivity can 
still be understood from Fig[T] Disorder leads to broadening 
of the energy levels and a finite density of states at the Dirac 
point. This implies that the Pauli exclusion is not effective in 
blocking transitions and hence there is always a finite conduc- 
tivity even form uj < 2p. The conductivity in the "forbidden" 
region increases with the increase in the number of impuri- 
ties. The fact that the imaginary part of the electron-phonon 
self-energy vanishes for electron energies between p — u>o and 
p + ujq indicates that for 2 x \p — ujq\ < ui < 2 x (ujq + p) 
the electron-phonon coupling does not produce any effect in 
the conductivity. For p < ujq, we expect the appearance of an 
anti-Stokes line at u>a.s. = 2 x (u>o + p) and a Stokes line at 
c^s. = 2 x (u>o — p)- For ujq < p the Stokes line lies inside of 
the Pauli blocked region and hence it should be suppressed. 

In Figure [3] we plot the infrared conductivity of a graphene 
in units of the far-infrared conductivity ctq = ire 2 /(2h). The 
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nf = 0.00013, n = 0.00004, T = 45 K 
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FIG. 3: (color online) Real part of the infrared conductivity including 
the effect of phonons, unitary scatterers and charged impurities. The 
parameters are T = 45 K, m = 4.0 x lCT 5 , and nf = 1.3 x 1(T 4 . 
The dashed vertical lines correspond to uj = 2/i and the shorter 
dotted-dashed to uj = 2(cjq + n) for different values of /i. 



main feature is that the conductivity is finite in the range 
< Tiw < 2/i and increases as the gate voltage decreases. 
We choose the concentration of unitary scatterers states in 
Fig. [3] to be one order of magnitude smaller than the one of 
Coulomb scatterers 16 , and therefore the conductivity is mainly 
controlled by phonons and charged impurities. Another fea- 
ture of the curves in Fig. [5] is the large broadening of the 
inter-band transition edge at fku = 2/i (indicated by vertical 
dashed lines). Note that this broadening is not due to tempera- 
ture but to charged impurities, instead. In fact, the broadening 
for all values of V g is larger when the conductivity is con- 
trolled by charged impurities. As expected, the coupling to 
phonons produces a anti-Stokes line centered at 2(uiq + /i). 
For gate voltages with /j, < u)q there appears a Stokes line 
at 2(u>q — fi). We find, however, that the Stokes line is very 
sensitive to disorder and is fast suppressed by the inclusion of 



charge impurities. The optical phonons thus induce a conduc- 
tivity larger than a around these frequencies. This effect is 
washed out at high temperatures and low frequencies. We also 
note that for large biases the conductivity in the Pauli-blocked 
region becomes weakly voltage dependent. All these effects 
seem to be consistent with the recent infrared measurements 
of graphene on a SiC<2 substrate-^. 

In this paper we have studied the infrared conductivity of 
graphene at finite chemical potential, generalizing the results 
of Ref. Jgj]. The calculation includes both the effect of dis- 
order (unitary scatterers and charged impurities) and the ef- 
fect of phonons. The effect of acoustic phonons is negligi- 
ble, since it induces an imaginary part of the electronic self- 
energy that is much smaller than the imaginary part induced 
by either impurities or optical phonons. We find that optical 
phonons and charge impurities produce important modifica- 
tions in the infrared absorption leading to large conductivities 
in the Pauli-blocked energy region of ui < 2/i. The optical 
phonons also produce a conductivity larger than ao around 
2/i < uj ~ 2(wo ± /i) due to Stokes and anti-Stokes pro- 
cesses. It is interesting to note that for frequencies away from 
the Dirac point the imaginary part of the self-energy due to 
optical phonons is linear in frequency, a behavior similar to 
that due to electron-electron interactions in graphene^. The 
most important approximation in our calculations is associ- 
ated with the fact that we have neglected completely fiexural, 
modes since we assume that they are pinned by the substrate 
and hence have very high excitation energy, that is, away from 
the infrared regime. We stress that the only free parameter in 
the calculation is the density of charge impurities, nf that can 
change from sample to sample. 
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